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Abstract: A new method for calibrating the hadron response of a segmented calorimeter is de- 
veloped and successfully applied to beam test data. It is based on a principal component analysis 
of energy deposits in the calorimeter layers, exploiting longitudinal shower development informa- 
tion to improve the measured energy resolution. Corrections for invisible hadronic energy and 
energy lost in dead material in front of and between the calorimeters of the ATLAS experiment 
were calculated with simulated Geant4 Monte Carlo events and used to reconstruct the energy of 
pions impinging on the calorimeters during the 2004 Barrel Combined Beam Test at the CERN H8 
area. For pion beams with energies between 20 GeV and 180 GeV, the particle energy is recon- 
structed within 3% and the energy resolution is improved by between 11% and 25% compared to 
the resolution at the electromagnetic scale. 
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1 Introduction 

In the general case of non-compensating calorimeters, the response to hadrons will be lower than 
the response to particles which only interact electromagnetically, such as electrons and photons. 
This is due to energy lost in hadronic showers in forms not measurable as an ionization signal, i.e., 
nuclear break-up, spallation and excitation, energy deposits arriving out of the sensitive time win- 
dow (such as delayed photons), soft neutrons, and particles escaping the detector [1-3]. Moreover, 
the calorimeter response will be non-linear, since a hadronic shower has both an electromagnetic 
and a hadronic component, with the size of the former increasing with shower energy [4]. In addi- 
tion, the large phase space of hadronic interactions leads to substantial fluctuations in the size of the 
electromagnetic shower component from event to event, degrading the measured energy resolution. 

ATLAS [5] is one of the multi-purpose physics experiments at the CERN Large Hadron Col- 
lider (LHC) [6]. Scientific goals include searching for the Higgs boson and looking for phenomena 
beyond the Standard Model of particle physics, such as supersymmetry. Many measurements to be 
performed by the LHC experiments rely on a correct and accurate energy reconstruction of hadronic 
final-state particles. In the central barrel region, the ATLAS calorimeters consist of the lead-liquid 
ai^gon (LAr) electromagnetic calorimeter and the Tile steel-scintillator hadronic calorimeter. Both 
calorimeters are intrinsically non-compensating. 

Various techniques for equalizing the electromagnetic and hadronic shower response, i.e., 
achieving compensation, have been proposed. For a review, see [3], chapter 3. Software-based 
offline calibration techniques can use the topology of the visible deposited energy to exploit spatial 
event-by-event information on shower fluctuations and derive energy corrections aimed at restor- 
ing linearity in the response and improving the energy resolution. For example, the calorimeter cell 
energy density has been used for the calorimeter in the HI experiment [7] and is planned to be used 
in ATLAS [8]. 

In this study, a calibration technique based on Monte Carlo simulation is developed to deal 
with compensating the response of a segmented calorimeter to hadrons and correcting for energy 
lost in the dead material between two calorimeter systems. The correlations between longitudinal 
energy deposits of the shower have been shown [9] to contain information on the electromagnetic 
and hadronic nature of the shower. This information is utilized by making a principal component 
analysis of the energies deposited in the different calorimeter layers. The calibration is applied to 
pion beam test data, taken at the 2004 ATLAS BaiTcl Combined Beam Test [10-14]. The method 
presented here is an alternative to the standard ATLAS calibration schemes. The application is quite 
specific to ATLAS, but the framework is general and it can be tested on any segmented calorimeter. 
Energy corrections based on the longitudinal shower development have been proposed by ATLAS 
in the context of jet calibration [15-17]. 

The following section explains the basic principles of the method. Section 3 details the ATLAS 
Barrel Combined Beam Test, while sections 4 and 5 discuss calibration to the electromagnetic scale 
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and event selection, respectively. The Geant4 Monte Carlo simulation used is described in section 
6. Then, section 7 gives the details of the implementation of the calibration method. In section 8, 
the method is validated based on Monte Carlo simulations of pions. In the Monte Carlo simulation, 
the effect of the compensation weights and the dead material corrections are evaluated separately. 
Lastly, the linearity and resolution of the final calibrated energy is considered. Section 9 discusses 
systematic uncertainties. Results of applying the method to real beam test data are presented in 
section 10. Finally, conclusions are drawn in section 11. 



2 The Layer Correlation method 

The Layer Correlation calibration method (LC in the following) is aimed at calibrating the response 
of a non-compensating longitudinally segmented calorimeter to hadrons. Exploiting the properties 
of hadronic showers to characterize fluctuations in the deposited invisible energy, it uses a prin- 
cipal component analysis [18] of the energy deposited in the calorimeter layers. Observables that 
describe the shower fluctuations should be able to discriminate between different connections to be 
applied to recover invisible losses due to hadronic interactions. Through the principal component 
analysis, it is possible to reduce the number of dimensions that the corrections depend on, while 
still capturing a large amount of event fluctuation information and maintaining a good separation 
between events with different content of invisible energy. 

To derive the connections, the interaction of the shower particles with the detector material 
is simulated with the Geant4 [19, 20] Monte Carlo simulation tooMt. In the simulation the true 
energy deposited in the calorimeters and the non-instrumented material is known. The covariance 
matrix between the calorimeter layer energy deposits is calculated. Diagonalizing it, a new orthog- 
onal basis in the space of layer energy deposits is derived. It consists of the eigenvectors of the 
covariance matrix. By sorting the eigenvectors in descending eigenvalue order, the projection of 
the energy deposits in the calorimeter layers along the first few eigenvectors are made to describe 
the most important fluctuations in the longitudinal shower development. 

Using this information, compensation weights — conecting for the non-compensation of the 
calorimeters — are derived in the form of two-dimensional look-up tables in the projections along 
the first two eigenvectors of the covariance matrix. One table is used for each calorimeter layer. 
The tables are thus functions of two different Unear combinations of the observed energy deposits 
in the layers. 

In addition, energy losses in non-instrumented material (so-called "dead material") will vary 
depending on the shower development. In the ATLAS banel region, these losses are primarily 
in the region between the LAr and Tile calorimeters. The eigenvectors of the covariance matrix 
considered above can also be used to conect for this, resulting in a unified treatment for compen- 
sation and dead material correction by deriving both corrections from the same set of observables. 
In this implementation, the dead material conections have an inherent dependence on the beam 
energy. This dependence is removed by employing an iteration scheme, where at each step the esti- 
mated energy of the former step is used, until the returned value is stable. A detailed mathematical 
description of the method is given in section 7. 
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Figure 1. The layout of the 2004 Combined Beam Test. 
3 The 2004 ATLAS Barrel Combined Beam Test 

The energy calibration procedure is applied to data gathered in the fall of 2004 during the ATLAS 
Barrel Combined Beam Test at the H8 beam line of the CERN SPS accelerator. A full slice of the 
ATLAS barrel region was installed (see Figure 1). This included, firstly, the inner tracker with the 
pixel detector, the silicon strip semiconductor tracker (SCT), and the straw tube transition radiation 
tracker (TRT); secondly, the LAr and Tile calorimeters; and thirdly, the muon spectrometer. The 
pixel and SCT detectors were surrounded by a magnet capable of producing a field of 2 T, although 
no magnetic field was applied in the runs used for this study. 

The pixel detector [5] comprises six modules, each consisting of a single silicon wafer with 
an array of 40 x 400 |am^ pixels. The modules were arranged in locations mimicking the ATLAS 
configuration, with an approximate angle of 20 degrees with respect to the incoming beam. The 
semiconductor tracker (SCT) [5] uses sets of stereo strips for tracking. Each module gives two hits, 
one in each direction. Eight modules, corresponding to those in the ATLAS end-cap, were used. 
The TRT [5] forms the outermost tracking system in ATLAS. It consists of a collection of 4 mm 
diameter polyimide straw tubes filled with a mixture of xenon, carbon dioxide, and oxygen [5]. 
Transition radiation is emitted when a charged particle crosses the interface between two media 
having different refractive index. The amount of emitted radiation depends on the Lorentz 7 factor 
of the particle. This makes it possible to discriminate between electrons and hadrons, given the 
much higher 7 factor of the former at a given energy, due to their smaller mass. 

Details of the ATLAS LAr electromagnetic calorimeter are described elsewhere [5, 21]. In 
the beam test one calorimeter module was used. The calorimeter is made from 2.21 mm thick 
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Figure 2. The layout of the 2004 Combined Beam Test. 



accordion-shaped lead absorbers glued between stainless steel cathodes. Three-layered anode elec- 
trodes are interleaved between the absorbers, spaced by 2 mm gaps over which a high voltage of 
2 kV is applied. The module was placed in a cryostat containing liquid argon. The signal is read 
out by capacitive coupling between the two outermost and the central layer of the anodes. In front 
of this accordion module a thin presampler module was mounted. It consists of two straight sectors 
with alternating cathode and anode electrodes glued between plates made of a fiber-glass epoxy 
composite (FR4). The Tile hadronic calorimeter consists of iron absorbers sandwiched between 
organic scintillator tiles. It is described in detail elsewhere [5, 22]. The tiles and absorbers are 
oriented parallel to the direction of incoming particles. Every cell of the calorimeter is read out by 
two wavelength-shifting fibers, which in turn aie grouped together and read out by photo-multiplier 
tubes (PMTs). 

The calorimeters were placed so that the beam impact angle corresponded to a pseudo-rapidity^ 
of T] = 0.45 in the ATLAS detector. At this angle, the expected amount of material in front of the 
calorimeters was about 0.44 A/, where A/ is the nuclear interaction length [3, 23]. This includes the 
LAr presampler. The LAr calorimeter proper is longitudinally segmented in three layers that extend 
in total for 1.35 A/. The dead material between the LAr and Tile calorimeters spans about 0.63 A/. 
Finally the three longitudinal segments of the Tile calorimeter stretched in total for about 8.18 A/. 
A sketch of this setup is shown in Figure 2. In total there are seven longitudinal calorimeter layers 
(the LAr presampler; the front, middle, and back layers of the LAr calorimeter; and the so-called 
A, BC, and D layers of the Tile calorimeter). The length of the individual calorimeter layers was 
0.32, 0.96, and 0.07 A/ in the LAr calorimeter and 1.61, 4.53, and 2.04 A/ in the Tile calorimeter. 

In addition, special beam-line detectors were installed to monitor the beam position and re- 
ject background events. Those include beam chambers monitoring the beam position and trigger 
scintillators. Beams consisting of electrons, photons, pions, protons, and muons were studied. In 
this analysis, pion beams with nominal momenta of 20, 50, 100, and 180 GeV were used (see Ta- 
ble 1). Data belong to the fully combined run period, where all detector sub-systems were present 

'atlas has a coordinate system centered on the interaction point, with the x axis pointing towards the center 
of the LHC ring, the y axis pointing straight up, and the z axis parallel to the beam. Pseudo-rapidity is defined as 
— ln(tan(0/2)), where 6 is the angle to the positive z axis. 
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and operational. No magnetic field was applied around the pixel and silicon strip detectors. The 
beams were produced by letting 400 GeV protons from the SPS accelerator impinge on a beryllium 
target, from which secondary pions are selected. For the run at 180 GeV, positrons were nominally 
selected after the target. However, the beam still contained a contamination of positively charged 
pions, which were selected and used for this analysis with the methods described in section 5.1. 



4 Calorimeter calibration to the electromagnetic scale 

4.1 Cell energy reconstruction 

The individual cells of the calorimeter are calibrated to the electromagnetic scale, i.e., with the 
aim of coiTcctly measuring the energy deposited in the cell by a purely electromagnetic shower. 
The calibration of electronics of the LAr calorimeter is described in detail in [24]. The method of 
optimal filtering [25] is used to reconstruct the amplitude of the shaped signal, which is sampled by 
an ADC (analog-to-digital converter) at 40 MHz. The amplitude is calculated as weighted sum of 
the samples, after a pedestal level measured using random triggers is subtracted. fpA^Mev//samp> 
a constant factor, converts the measured current to an energy measured in MeV. The energy de- 
posited in the lead absorbers is taken into account by the sampling fraction /samp- The shaping 
electronics is calibrated by inserting calibration pulses of known amplitude. In the Tile calorimeter 
a parameterized pulse shape is fitted to the samples. A charge injection system is used to calibrate 
the read-out electronics, while a cesium source is used to equalize the cell response, including the 
response of the PMTs (see, for example [26]). 

4.2 Topological clustering 

Calorimeter cells calibrated to the electromagnetic scale are combined by adding up the energy in 
neighboring cells using a topological cluster algorithm [27]. The algorithm has thi^ee adjustable 
thresholds: Seed (S), Neighbor (A^), and Boundary (B). First, seed cells having an energy above the 

5 threshold are found and a cluster is formed starting with this cell. Then, neighboring cells having 
an energy above the threshold are added to the cluster. This process is repeated until the cluster 
has no neighbors with an energy above the A'^ thr^eshold. Finally, all neighboring cells having an 
energy above the B threshold are added to the cluster. To avoid bias, the absolute values of the cell 
energies are used. The S, N, and B thresholds are set to, respectively, four, two, and zero times the 
expected noise standard deviation in the cell considered. 

4.3 Pion energy reconstruction 

The reconstructed energy in a calorimeter layer L is obtained by considering all the topological 
clusters in the event and summing up the parts of the clusters that are part of that calorimeter layer. 
The total reconstructed energy is then derived by summing over the A'^iay longitudinal layers in the 
calorimeter. 
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^b"eam(GeV) 


£'meas (GcV) 


No. ev. bef. cuts 


No. ev. after cuts 


/prot 


20 


20.16 


49871 


8957 


< 17% (84% CL) 


50 


50.29 


109198 


29578 


( 45 ± 12)% 


100 


99.89 


67220 


5843 


( 61 ± 6)% 


180 


179.68 


105082 


11780 


( 76 ± 4)% 



Table 1. Data samples taken in the 2004 Combined Beam Test used in the present analysis. 

5 Event selection and particle identification 

5.1 Event selection 

A signal in the trigger scintillator and a measurement in adjacent beam chambers that is compatible 
with one particle passing close to the nominal beam line are required. In addition, exactly one track, 
where the sum of the number of hits in the Pixel detector and the SCT is more than six, is asked 
for, as well as at least 20 hits in the TRT. The track in the TRT must be compatible with a pion 
track, i.e., no more than two hits passing the high threshold must be present. Events with a second 
track in the TRT are rejected: this ensures that the pion does not interact strongly before the TRT. 
Furthermore, there must be at least one topological cluster (see section 4.2) with at least 5 GeV 
in the calorimeter. This cut rejects muons contained in the beam and does not influence the pion 
energy measurement. To reject some residual electron background, events with more than 99% of 
their energy in the LAr calorimeter are excluded. The same selection is applied on simulated Monte 
Carlo events as on data, with the exception of cuts related to the beam chambers and scintillators. 

5.2 Proton contamination 

This study used beams of pions with positive electric charge. These beams are known to have a 
sizable proton contamination /prot defined as the fraction of events in a sample that result from 
protons impinging on the calorimeters. It varies between different beam energies. The TRT makes 
it possible to measure the average proton contamination of the test beam for each beam energy, 
owing to the different probabilities between pions and protons of emitting transition radiation, 
although it is not possible to discriminate between the particles on an event-by-event basis. The 
measured [10] contamination is reported in Table 1. For the 20 GeV beam energy, a one-sided 
confidence interval is given. In the analysis, a proton contamination of 0% was used. Agreement 
is found with measurements performed by a Cerenkov counter at a 2002 beam test [28] conducted 
in the same beam line. 

6 Monte Carlo simulation 
6.1 Hadronic shower simulation 

All calibration corrections are extracted from a Geant4.7 [19, 20] Monte Carlo simulation, with an 
accurate description of the Combined Beam Test geometry. The physics list — i.e., set of models — 
QGSP_BERT was used. It uses the QGSP [29] (Quark Gluon String Pre-compound) phenomeno- 
logical model describing the hadron-nucleus interaction by the formation and fragmentation of 
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excited strings together witii tiie de-excitation of an excited nucleus. The Bertini model [30-32] 
of the intra-nuclear hadronic cascade is used to describe nuclear interactions at low energies. This 
model treats the particles in the cascade as classical and propagates them through the nucleus, 
which is modeled as a medium with a density averaged in concentric spheres. Excited states aie 
collected and the nucleus decays in a slower phase following the fast intra-nuclear cascade. 

The Bertini model is applied up to an energy of 9.9 GeV, while the QGSP model applies from 
12 GeV and upward. In an intermediate range of 9.5-25 GeV, the low-energy parameterized LEP 
model [33] is used. In the energy ranges where models overlap, the decision which one to use is 
made stochastically using a continuous linear probability distribution that goes from exclusively 
using the low-energy model at the lower end of the region to exclusively using the high-energy 
model at the upper end. 

6.2 Detector simulation 

The simulation provides not only reconstructed calorimeter cell energies at the electromagnetic 
scale — including the effects of the readout electronics — but also the true deposited energy, which 
is divided into four components: electromagnetic visible, hadronic visible, invisible, and escaped. 
Visible energy results from ionization of the calorimeter material. Invisible energy is energy not 
directly measurable in the detector, such as break-up energy in nuclear interactions. The escaped 
energy represents the small contribution from neutrinos, high-energy muons and, possibly, neutrons 
and low-energy photons escaping the total simulated volume. 

6.3 Event samples 

Monte Carlo samples were produced by simulating both pions and protons impinging on the de- 
tector setup. Two statistically independent event samples were produced by dividing the available 
sample into two approximately equal parts: one set ("correction" samples in the following) was 
used to derive compensation weights and dead material corrections, while the other set ("signal" 
samples in the following) was used to validate the weighting procedure and find the expected per- 
formance. Pions and protons were simulated at 25 different beam energies, ranging from 15 GeV 
to 230 GeV. In total, about 800 000 events per sample and particle type were available after event 
selection. The energy spacing was 2, 3, or 5 GeV up to 70 GeV and 10 or 20 GeV above 70 GeV. 
This spacing was found to give satisfactory performance (see sections 8 and 10). Further studies of 
different spacings can be pursued when applying this technique to different calorimeters to explore 
possible improvement in performance. 

Taking the proton beam contamination mentioned in section 5.2 into account, all the available 
"correction" Monte Carlo samples were used to build a "mixed" pion-proton sample, one for each 
energy available in the data (see Table 1). Each of these samples is used as input when deriving 
the corrections used for that proton fraction. In this way the corrections were tuned to the studied 
proton fraction. If the samples had different numbers of events, a sample-dependent weight was 
first applied to give them equal weight before selection cuts. Then, given the proton contamination 
/prot at a given energy, pion and proton events for each same-energy pair of samples were assigned 
a weight of 1 — /prot and /prot, respectively. 
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7 Implementation of the Layer Correlation method 
7.1 Calculation of the eigenvectors of the covariance matrix 

Each event is associated with a set of A'^iay layer energy deposits (E\^'^, ... , one per calorimeter 

layer, representing a point in an A'^iay -dimensional vector space, referred to in the following as 
the space of layer energy deposits. They are reconstructed energies at the electromagnetic scale, 
formed as calorimeter layer sums of topological clusters as described in section 4.1. The May- 
dimensional covariance matrix of the layer energy deposits is calculated as 

Cov(M,L) = {E^^^ED - {ED , (7.1) 

where M and L denote calorimeter layers and £'^'^ is the energy reconstructed at the electromagnetic 
scale in calorimeter layer M. The averages are defined as 

{E'^'ED =^=^^ and {E^ =^^. (7-2) 

The sums are performed over all the Mv events in the sample. The eigenvectors of the covariance 
matrix form a new orthogonal basis in the space of layer energy deposits. The coordinates of the 
point in the May -dimensional vector space coiTesponding to an event / can be expressed in this new 
eigenvector basis as 

■^eig,M = 12 (^M^.lEl^, (7.3) 
L 

where (x}^\ are the coefficients of the transition matrix to the new basis. Projections of events 
along the covariance matrix eigenvectors represent independent fluctuations. The variances of 
those fluctuations are given by the corresponding eigenvalues. The eigenvectors are sorted in de- 
scending order according to their eigenvalues, meaning that the first eigenvectors determine the 
directions along which most of the event fluctuations take place. The layer energy covariance 
matrix Cov(M,L) (equations 7.1 and 7.2) is calculated using events from the "mixed" sample. 

In any given event a symmetric energy cut is applied on each layer energy such that the energy 
for that layer is re-defined as E^j^^, if \E^^^\ ) Ef^', zero otherwise. The goal of such cuts is to elimi- 
nate the contribution of noise-dominated layers. The energy threshold values for each calorimeter 
layer can be found in Table 2. The cuts were optimized to obtain the best expected compensation 
performance on Monte Carlo samples at 50 GeV. 

A physical interpretation of the eigenvalues and normalized eigenvectors can be obtained from 
Figure 3, which shows the components of the first three eigenvectors expressed in the original basis 



of calorimeter layer energy deposits. We find that 

■^efg.O ~ -^(— 2£'LAr,middle +£"1116^ + -^Tile,Bc), (7.4) 

KHi - -^i-EjacA + £Tiie,Bc) , and (7.5) 

■^efg,2 ~ (£'LAr,middle + £"1116^ + ^'xile.Bc) • (7-6) 
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Calorimeter layer 


I nresnoiQ (yjc v ) 





0.032 


1 


A 1 AO 

O.lOo 


2 


0.030 


3 


0.150 


4 


0.039 


5 


0.070 


6 


0.042 



Table 2. Energy thresholds per calorimeter layer 




Calorimeter samplings 

Figure 3. Eigenvector components for the first three eigenvectors expressed in the basis of the seven layers 
of the ATLAS calorimeters in the Combined Beam Test for a simulated mix of protons and pions with 45% 
proton contamination. 



So in a qualitative but suggestive way, we can make the interpretation that f'gigo corresponds to 
the difference between the Tile and LAr calorimeters, since most of the energy deposited in the 
LAr calorimeter is deposited in the middle layer. Elf^ j corresponds to the difference between the 
second and first layers of the Tile calorimeter, while E^f^ 2 corresponds to most of the energy of the 
event. The other eigenvectors represent individual calorimeter layers. These layers are rather thin 
and appear to be uncorrelated to the other layers. 
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7.2 Compensation weights 

The compensation weights account for the non-linear response of the calorimeters to hadrons. 
There is one weight table for each calorimeter layer, i.e., three for the LAr calorimeter and three 
for the Tile calorimeter. The seventh layer, the LAr presampler, which in order is the first layer, is 
not used in the weighting procedure, as explained below. The total reconstructed energy is the sum 
of the weighted energies in each calorimeter layer: 

^weighted ^ ^^^rec (7 7^ 

^weighted _ ^ ^weighted ^ 
L 

For each event /, there is an ideal set of A'^iay coefficients that would re-weight each recon- 
structed energy deposit in layer L to the true deposited energy: 

ideal _ ptrue //rrec n\ 

The symbol Ej^f (Ef^^) denotes the reconstructed (true) energy deposited in the layer in the fi^ 
event. The task is to find a set of weights wl that approximate the ideal weights. In general, for 
each layer L, the weight is an A'lay -dimensional function of the layer energy deposits. Exploiting 
the fluctuation-capturing properties of the eigenvector projections, the weights can in general be 
derived as a function of an A/^-dimensional subspace of the A/^iay-dimensional space of layer energy 
deposits, spanned by the first A'^ eigenvectors. In the absence of an analytic formulation, the layer 
weights Wl are estimated by Monte Carlo sampling: multi-dimensional cells are built, which parti- 
tion the A^-dimensional vector space along the directions of the base eigenvectors. In general, these 
cells are multi-dimensional hyper-cubes. They are referred to as bins below. 

For each bin k one defines the weight as the average of the ideal weights of equation 7.9: 

w,,L = {E'Ci/El'Dk = 77^1^07^1!^ (7-10) 

where the summation is performed for the A^evi events in the bin. If each event has a weight^ pi, 
the average is modified accordingly: 

y. r,.£'true /^rec 

w,,, = {E^VEl^JU = (7.11) 

Li Pi 

Using bin k of the weight tables, the total reconstructed energy becomes 

Ez:r'=j:w,,LEr. had 

L 

Here, the wic.l functions defined in equation 7.1 1 are estimated in bins of the two-dimensional space 
spanned by the eigenvectors corresponding to the two highest eigenvalues, i.e., N = 2. Thus each 
layer is associated with a two-dimensional look-up table. For a given layer the average weights in 
each two-dimensional bin are calculated using only the energy values that passed the cuts defined 
in section 7.1. The table has the same number of equally spaced bins along the two dimensions: 



For instance, to equalize the number of events for all data sets. 
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128 X 128. Bi-linear interpolation is performed between the bins. Weights for the LAr presampler 
are not calculated, even if the presampler is kept in the covariance matrix. No weights are applied 
to the energy deposited in the presampler layer, and energy deposited in the presampler itself is 
taken as part of the upstream dead material losses. 

In addition the compensation weights and corrections derived from the proton sample are 
corrected by the factor 

pnom 
beam 



■^beam "'proton 



where wiproton is the proton mass, to account for the fact that, for a proton, the sum of the total true 
deposited energy in the calorimeter is E^°^ — wiproton- 

Typical compensation weight tables are shown in Figure 4: they illustrate the look-up tables 
for the second (middle) layer of the LAr calorimeter and for the first and second layer of the Tile 
calorimeter for a pion-proton mixed sample with 45% contamination. The triangular shape visible 
in the weight tables can be understood from the interpretation of the eigenvectors of equations 7.4 
and 7.5. With increasing energy in the Tile calorimeter and less in the LAr calorimeter, i.e., i'gf^o 
is large, there are more values that can be assumed by Elf^ , , which is the approximate difference 
between the first and second layers of the Tile calorimeter Three lines can be seen extending from 
the origin to each of the three corners of the triangle. Firstly, the line extending from the origin and 
to the left corresponds to events where close to all of the energy is deposited in the LAr calorimeter. 
The small slope is due to the slight dependence of E^?^ j on the second layer of the LAr calorimeter 
Secondly, the line extending up and to the right corresponds to events where all energy is deposited 
in the second layer of the Tile calorimeter. Along that line, weights aie small for the first sampling 
of the Tile calorimeter, since particles aie still minimum-ionizing in that layer. Thirdly the faint line 
extending down and to the right con^esponds to events where close to all the energy is deposited in 
the first layer of the Tile calorimeter 



7.3 Dead material corrections 

Regions of dead material constitute those parts of the experiment that are neither active calorimeter 
read-out material (liquid argon or scintillator), nor sampling calorimeter absorbers (mostly lead or 
steel). The LC technique is used for the dead material between the LAr and the Tile calorimeters, 
while a simple parameterized model is utilized for other losses. 

7.3.1 Dead material between tlie LAr and Tile calorimeters 

Most of the dead material is in the LAr cryostat wall between the LAr and Tile calorimeters. In 
this 0.6 A/ region, pion showers are often fully developed, giving rise to large energy loss. Each 
event / is associated with a point in the layer energy deposit vector space as explained in section 
7. 1. It also has a true total energy lost in the dead material between the LAr and Tile calorimeters: 
^LArTiiT(')- The dead material correction £"^1116 ^^'^ ^^^^ event / can be derived as a T-dimensional 
function of the layer energy deposits. In general, the subspace chosen for deriving the dead material 
correction and its dimension T can be different from the one chosen for compensation, both in 
content (spanned by different eigenvectors) and in dimension {T can be different from A'^). The 
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(a) 




(c) 



Figure 4. Compensation weights as a function of the first two eigenvector projections for simulated pion- 
proton mixed events (45% proton contamination) in the second layer of the LAr calorimeter (a), first layer 
of the Tile calorimeter (b), and second layer of the Tile calorimeter (c). 
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Figure 5. Look-up table for LAr-Tile dead material corrections as a function of the first and third eigenvector 
projections normalized to beam energy for 45% proton contamination. 

value of £^J^^xiie estimated by Monte Carlo sampling. For any T-dimensional bin m one defines 

pDM _ /pDM.true \ /-y i 

■^LArTile.m " X-^^LArTile,!/™ ' ^'■'^'*) 

where the average is performed for the events in that bin. 

Here, the correction defined in equation 7.14 is calculated in bins of the two-dimensional 
space spanned by the eigenvectors corresponding to the first and third eigenvalues, i.e., T = 2. 
This was the combination of eigenvectors that was found to give the best performance. As for 
the compensation weights, correction tables are derived from a 128 x 128 bin look-up table and 
bi-linear interpolation is performed between the bins. 

The three dimensions of the look-up table are all shown to scale with the beam energy, i.e., 
a table determined at a given beam energy can be turned into one at a different beam energy by 
scaling all the dimensions with the ratio of the two energies. Consequently, all dimensions in the 
table — the eigenvector projections and the average dead material losses — are divided by the beam 
energy when filling the table. That is, the event coordinates in the space of layer energy deposits 
are expressed as 

^d™ = EZm/E = j:al!l,Er/E, (7.15) 

L 

where the variables have the same meaning as in equation 7.3 and E is the best estimate of the 
beam energy of the simulated pion in that event (see below). The dead material look-up table is 
shown in Figure 5 for a pion-proton mixed sample with 45% contamination. The figure shows the 
distribution of the rescaled dead material energy as a function of the rescaled event coordinates. 
Regions with different dead material fractions can be differentiated. They range between and 
more than 30% of beam energy. In addition, the samples at different energies behave very similarly 
as a function of the re-scaled variables. 
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Figure 6. Mean dead material losses other than those between the LAr and Tile calorimeters as a function 
of the beam energy. Filled circles indicate the mean loss obtained from Monte Carlo simulation. The line 
indicates a parameterization to interpolate between the beam energies. 



7.3.2 Other dead material corrections 

While the energy losses between the LAr and Tile calorimeters dominate, there are still other 
regions where dead material losses can occur. These are losses located in the material upstream of 
the LAr calorimeter, between the LAr presampler and the first LAr calorimeter layer, and energy 
leakage beyond the Tile calorimeter. To compensate for these losses the mean energy loss was 
determined as a function of beam energy and the resulting data points were fitted using a suitable 
functional form 

cDM / 77 \ J C[+C2 \/ £'beam if £'beam < Eq 

^other I-beam j — i ^ , ^ /c- 17 \ ■ C'-lOj 

[C3 + C4 (£beam " Eo) Otherwise, 

where Eq = 30 GeV. As an example, the fit for a proton fraction of 45% can be seen in Figure 6. 
The resulting fitted parameters are 

Ci = (-353 ± 23) MeV, (7. 17) 

C2 = (8.47 ±0.17) VMcV, (7.18) 
C3 = (1102±3) MeV, and (7.19) 
C4 = 0.01392 ±0.0001. (7.20) 
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7.4 Applying the calibration 

The final energy after calibration consists of the sum of the weighted calorimeter layer energies 
and the dead material coiTcctions: 

E-coiT _ weighted ,pDM(p^ .7 9 H 

The index k stands for the bin in the appropriate A'^-dimensional space of layer energy deposits used 
in the weight tables (equations 7. 10 or 7. 1 1), while m is the bin in the T-dimensional space of layer 
energy deposits used to build the LC estimate for the energy loss in the dead material between 
the LAr and Tile calorimeters obtained from equation 7.14. The total dead material correction is 
derived from summing the two contributions derived in sections 7.3.1 and 7.3.2: 

^tot.m \^ ) - -^LArTile.m + ^other J ' {1.11} 

where E is the best estimate for the total deposited pion energy used to estimate ^beam in equa- 
tion 7.16. 

The events in a Monte Carlo sample are usually generated at a fixed beam energy in order 
to test the calorimeter response. Corrections derived from a fixed beam energy sample are, in 
principle, dependent on that information, i.e., they depend on the same quantity (pion energy) for 
the reconstruction of which they should be used. For the compensation weights, this dependence is 
overcome by superposing events from all the available energies. The eigenvector projections scale 
approximately with the energy of the incoming particle, meaning that regions in the table that come 
in use for a certain particle energy will be dominated by samples close to that energy. 

On the other hand, the look-up-table-based LAr-Tile dead material correction and the param- 
eterized model for the other dead material losses have an inherent dependence on an assumed beam 
energy when applying the corrections (see equations 7.15 and 7.22). This dependence is overcome 
using an iteration technique, giving the end result of depending only on the energy in the calorime- 
ters. At each step the best estimate of the reconstructed energy ^toT^ after all corrections is used to 
set both the scaling factor \/E (equation 7.15) for the LAr-Tile conection and the best pion energy 
estimate in the parameterization for the other dead material corrections. Each new estimate of the 
energy is used to pick up a new correction from the look-up table until the returned value is stable. 
In the initial step E^°" is just the pion energy after compensation weights ai^e applied. The iteration 
cut-off is a tunable parameter. 

The process of applying the calibration is as follows: 

• Associate each event to a bin in both the A/^-dimensional compensation weight and the T- 
dimensional dead material correction spaces defined in sections 7.2 and 7.3 by expressing its 
electromagnetic-scale energy deposit vector in the new eigenvector basis derived from the 
simulated events. 

• Extract compensation coiTcctions for the energy of each given layer and the LAr-Tile dead 
material correction from the look-up tables. Apply all coiTcctions according to equations 7.21 
and 7.22. 

• Use the iteration for dead material corrections. 
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8 Method validation on Monte Carlo simulation 



Before applying it to beam test data, the calibration is validated on a Monte Carlo sample statis- 
tically independent of the one used for extracting the corrections. First, the performance of the 
compensation weights is evaluated, then the linearity and resolution of the method as a whole. The 
weighting technique is validated on Monte Carlo simulation samples in separate steps: 

• Reconstruct the true deposited energy in the calorimeters (compensation validation). 

• Reconstruct the full energy of the incoming particles, including dead material corrections, 
and quantify the performance in terms of hnearity and resolution. 

The performance is evaluated in terms of bias and resolution. The weights and dead material 
corrections are derived from the "correction samples" and applied on the statistically independent 
"signal samples" (see section 6.3). The results in this section are derived for pions only. 

8.1 Compensation validation 

The reconstructed pion energy after compensation correction is compared to the true deposited en- 
ergy in the calorimeter. The event-by-event difference e^^^^^^^^ — £''™'^(calo) is considered, where 
£'j™''(calo) is the true total energy deposited in the calorimeter. The bias in the energy recon- 
struction is defined as the average value {E^i^^^^^'^ — £'J™'^(calo)) and the resolution is obtained by 
calculating the standard deviation a{E^^^^^^^^ — £'J™'^(calo)). 

The performance of the LC technique is compared with a simple calibration scheme (called 
/comp in the following) which uses beam energy information: each event in the sample is weighted 
with the same factor /comp = (£'t™'')/(^toD' where (£'t™'') ((f'toD) is the average true total (re- 
constructed) energy deposited in the given sample in the whole calorimeter, but not in the dead 
material. The /comp calibration scheme provides a reference scale to which the improvement in 
resolution of the LC weighting can be compared. 

The results of the validation procedure are shown in Figure 7. By construction, there is no bias 
in the energy reconstruction for the calibration procedure using a simple factor. The LC weighting 
mostly gives a slight positive bias of about 0.6%. At the lower edge of the energy range studied, 
the bias instead turns slightly negative. The resolution improvement increases with beam energy. 
It is about 10% at 50 GeV and about 20% at 180 GeV. 

8.2 Dead material corrections 

Figure 8 shows the bias of the weighted energy, and also the bias of the dead material corrections. 
For most energies, the LAr-Tile dead material correction has a slight negative bias, while at low 
energies the bias is positive. The bias is 0.5% maximally. This cancels out most of the bias from 
the weighting. The final energy is reconstructed correctly within a few per mil. 

8.3 Linearity and resolution in the Monte Carlo sample 

The performance for the fully corrected energy reconstruction is finally assessed in terms of linear- 
ity with respect to the beam energy and relative resolution. The reconstructed energy distribution 
is fitted with a Gaussian distribution in the interval (/i - 2a, jJ. + 2a), where jj. and a are the mean 
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Figure 7. Bias (a) and resolution (b) of the reconstructed energy after compensation correction minus 
the true deposited energy for energy deposited in the calorimeters in simulated samples for the calibration 
procedure using a simple factor and LC weighting. 



value and the standard deviation, respectively. This interval is found iteratively. The mean value 
Efii and the standard deviation agt of the fitted Gaussian are used together with the beam energy 
£^beam to define the linearity and the relative resolution. 
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Figure 8. Bias (reconstructed energy minus true deposited energy, divided by beam energy) for the three 
individual corrections: Weighted calorimeter energy, correction for dead material between the LAr and Tile 
calorimeters, and other dead material corrections. Lastly, the bias of the final reconstructed energy, which is 
the sum of the three. 

• The linearity is £'fit/£'beam as a function of ^beam- 

• The relative resolution is afit/^fit as a function of ^beam- 

Both linearity and relative resolution are derived for the energy distribution at four stages of the 
energy reconstruction: 

• at the electromagnetic scale, 

• after applying the compensation weights, 

• after compensation weights and application of dead material correction for losses between 
the LAr and Tile calorimeters, and 

• after compensation weighting and all dead material corrections. This last step aims to recon- 
struct the pion energy. 

This is shown in Figure 9. At the electromagnetic scale the calorimeter response is non- 
linear — as expected — and only about two thirds of the pion energy is measured. After weighting, 
between 80% and 90% of the incoming pion energy is recovered, while the dead material between 
the LAr and Tile calorimeters accounts for an additional 8-10%. After all corrections the correct 
pion energy is reconstructed within 1% for all beam energies. Each correction step makes the 
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calorimeter response more linear. The compensation weights give a better improvement of the lin- 
earity at high energies, while the dead material effects play a more significant role at low energies, 
in particular at 20 GeV where other corrections than LAr-Tile dead material are important to get 
to within 1% of the beam energy. The relative resolution is improved when applying each of the 
different correction steps. ^ At high beam energies (above ^beam = 100 GeV) the contribution of the 
compensation weights to the improvement in energy resolution has the same magnitude as that of 
the LAr-Tile dead material corrections. At lower beam energies dead material coiTcctions account 
for about 70% of the relative resolution improvement down to about ^beam — 30 GeV. Below ^beam 
~ 30 GeV all the connections account for a similar fraction of the improvement: other dead material 
corrections than those for LAr-Tile account for about 20% of the resolution improvement, they are 
marginal above that threshold. 

9 Systematic uncertainties 

Systematic uncertainties on the calibrated energy can be divided into 

• The uncertainty of the beam energy: 0.7% [13]. 

• The absolute electromagnetic scale uncertainty, which is estimated to be 0.7% [13] in the 
LAr calorimeter and 1.0% [26] in the Tile calorimeter. Scaling the cell energies with their 
corresponding uncertainties gives a combined electromagnetic scale uncertainty of 0.9%. 

• The sensitivity of the results to the proton fraction at each beam energy. It was estimated by 
varying the fraction used to calculate the connections. With the assumed fraction adjusted up 
or down one standard deviation of the TRT measurement, the relative variation in linearity 
and resolution in data and Monte Carlo simulation was found to be of the order of 1% for 
£'beam = 20 GcV and 50 GeV and less than 0.5% for 100 GeV and 180 GeV. 

Adding these contributions in quadrature gives a total systematic uncertainty of less than 2% for 
each beam energy. 

10 Application of the method to beam test data 

Finally, the method is applied to beam test data, which is compared with Monte Carlo samples with 
a weighted mixture of pions and protons to match the beam composition. 

10.1 Data to Monte Carlo simulation comparison 

The pion-proton "mixed signal" samples are used to compare data and Monte Carlo simulations 
in terms of the distribution of the first three components of the layer energy vector along the basis 
of covariance matrix eigenvectors as defined in section 7.1. Figure 10 shows such a comparison 
for a proton fraction of 45% and a beam energy of 50 GeV. Good agreement is obtained between 
data and simulation. The distribution for ^eig o shows a double peak structure that separates events 
mainly showering in the Tile calorimeter from those where the shower starts earlier. 

^The apparent discontinuity in resolution between the results at energies below 150 GeV and those above might be 
due to a geometry change in the description of the beam test setup: three centimeters of aluminum were included in the 
Inner Detector system for energies larger than or equal to 150 GeV. 
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Figure 9. Linearity (a) and relative resolution (b) for simulated pure pion samples at the electromagnetic 
scale, with compensation weights applied, compensation weights plus LAr-Tile dead material correction 
applied, and all corrections applied. 

The shapes of the energy distributions (in unit bins of energy and events) for data and Monte 
Carlo simulation are compared in Figure 1 1 . The corrections are successively applied. Already 
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Figure 10. Distribution of the first three eigenvector components for data (filled circles) and Monte Carlo 
simulation pion-proton "mixed signal" with a proton fraction of 45% and a beam energy of 50 GeV. 

at the electromagnetic scale the energy distribution is not well reproduced. The distribution in the 
Monte Carlo simulation is narrower and less skewed than in the data. This effect is even larger at 
20 GeV but less pronounced at higher energies. The quality of the initial description of data by 
Monte Carlo simulation is not modified by the application of the compensation weights and dead 
material corrections (see also section 10.2). 

10.2 Linearity and resolution on data 

The performance of the method, as applied to simulation and real beam test data is shown in Figure 
12. The data at the electromagnetic scale in this analysis and the one presented in references [11, 
12] are in reasonable agreement. The largest deviations of about 3% (2%) are seen at 20 GeV 
(180 GeV). At 20 GeV the difference can be explained by the fact that in that study, the energy in 
the calorimeters — instead of using topological clustering — was determined by adding the energies 
^ceii of those calorimeter cells having a pseudo-rapidity within 0.2 of the beam impact point and for 
which ^ceii is two standard deviations above the expected noise. At 1 80 GeV, data in that study were 
taken with a beam of negatively charged pions, which does not suffer from proton contamination. 
In addition, for all beam energies data were taken in an earlier run period with a different material 
configuration upstream of the calorimeters. 

The linearity and relative resolution are extracted at all beam energies for both data and "mixed 
signal" Monte Carlo samples. As in section 8.3 the reconstructed energy distribution is fitted with 
a Gaussian distribution in the interval (jj. - 2a, jJ. + 2a), where jj. and a are the mean value and 
the standard deviation, respectively. Data (simulation) are shown with markers (horizontal lines) 
at the electromagnetic scale, with compensation weights applied, with the dead material correction 



-22- 



• EM scale (data) 
■ Weighted (data) 

A Weighted + LAr-Tiie DM corr. (data) 



— EM scale (MC) 

— Weighted (MC) 

Weighted + LAr-Tile DM corr. (MC) 



T Weighted + all DM corrs. (data) — Weighted + all DM corrs. (MC) 



> 
CD 

cs 



LJJ 
T3 



T 



T 



T 



T 



T 



T 



0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 





60 70 

E (GeV) 



Figure 11. Normalized energy distribution for fibeam = 50 GeV after applying subsequent corrections for 
compensation and dead material effects. For the Monte Carlo simulation a proton fraction of 45% is used. . 



for energy lost in dead material between the calorimeters applied, and lastly at the final calibrated 
stage, including all dead material corrections. 

After all calibration connections, the linearity is recovered within 3% for all beam energies. 
The discrepancies between data and Monte Carlo are inherited from the reconstructed energy at the 
electromagnetic scale and they are not considerably changed when the calibration is applied. The 
relative resolution in data is improved by about 1 1% at low energy (20 GeV) and about 25% at high 
energy ( 1 80 GeV) when moving from the electromagnetic scale to the fully corrected energy scale. 
A similar relative improvement is obtained in the Monte Carlo simulation: 14% at low energy and 
24% at high energy. The relative resolution is, however, smaller in Monte Carlo simulation than 
in data: the discrepancies, at each correction stage, vary between 10% and 21% depending on the 
energy. The relative resolution is smaller in Monte Carlo simulation than in data already at the 
electromagnetic scale, by about 10-16%, depending on beam energy. The discrepancies in the 
shape of the total energy distribution are more pronounced at lower energies and they are already 
present at the electromagnetic scale. 

The ratio of data to Monte Carlo simulation is unchanged within 1% (4%) for linearity (res- 
olution) after the corrections are apphed. For linearity such changes are of the same order of 
magnitude as the discrepancies between data and Monte Carlo simulation at the electromagnetic 
scale: the agreement between data and Monte Carlo simulation is the same for all correction stages. 
This means that the Monte Carlo simulation is able to predict the corrections that should be applied 
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Figure 12. Data and Monte Carlo simulation are compared for linearity (a) and relative resolution (b) at all 
stages of the corrections. The markers show data, while the horizontal lines show (a) from bottom to top and 
(b) from top to bottom: electromagnetic scale, compensation weights applied, weight plus LAr-Tile dead 
material correction applied, and aU corrections applied. See text for details. 
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on the data. 

The abiUty of the Monte Carlo simulation to reproduce the data at the electromagnetic scale 
(i.e., before any coiTcction) seems to be the most critical limiting factor. For the relative resolu- 
tion the changes are small, if compared with the discrepancies at the electromagnetic scale: the 
discrepancies do not get worse when the corrections ai^e applied to the data. From preliminary 
studies a newer Geant4 version (4.9) is able to provide a better, but still not good, description of 
the resolution in the data. 

11 Conclusions 

An energy calibration technique was developed to deal in a coherent manner with both compensat- 
ing the hadron response and correcting for the most significant dead material losses in a segmented 
calorimeter. The technique is based on the sensitivity of the correlation between the deposited 
energies in the different calorimeter layers to hadronic and electromagnetic deposits. 

The calibration technique was successfully applied to the energy reconstruction of pions im- 
pinging on a subset of the central ATLAS calorimeters during the ATLAS combined beam test in 
2004. When taking into account the beam composition of pions and protons, linearity is recov- 
ered within 3% and relative resolution is improved by between 11% and 25%. Consistency with 
the expectation from Monte Carlo simulation studies is good for both the linearity and the per- 
centage improvement in relative resolution. The absolute value of the relative resolution (after all 
corrections) is larger in data than Monte Carlo simulation by 10% to 21%. 

The discrepancies between data and Monte Carlo simulation are inherited from the recon- 
structed energy at the electromagnetic scale and they are not considerably altered when applying 
the calibration. Additional improvement in the data description by Monte Carlo simulation can 
help to fulfill the expected absolute value for the relative resolution. 
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